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Abstract 

We consider supersymmetric models in which the neutral Wino is the lightest su- 
perparticle (LSP), and study the possibility that non-thermally produced Wino plays 
the role of dark matter. The thermal relic density of Wino is smaller than the present 
mass density of dark matter if < 2.9 TeV; however, even with smaller Wino mass, 
Wino can be the dominant component of dark matter if it is non-thermally produced 
by the decay of an exotic particle. In order to study such a possibility in detail, we 
perform a precise calculation of the present mass density of Wino produced by the 
exotic-particle decay, taking account of the following effects: (i) coannihilation among 
charged and neutral Winos, and (ii) Sommerfeld effect on the pair annihilation cross 
section of Winos. We consider several well-motivated cases where the exotic particle 
corresponds to cosmological moduli fields, gravitino, or axino, and discuss the impli- 
cation of the Wino LSP for these cases. 



1 Introduction 



In models with low-energy supersymmetry (SUSY), there may exist various "exotic" particles 
(which we call X) which have very weak interactions with the particles in the minimal SUSY 
standard model (MSSM). The examples include moduli fields, Polonyi field (responsible for 
SUSY breaking), gravitino (which is the superpartner of graviton), axino (which is the 
superpartner of axion in SUSY Peccei-Quinn model), and so on. In most of the cases, 
these particles are irrelevant for low-energy phenomenology because of the weakness of their 
interactions. 

Cosmologically, however, they often cause serious problems pQ |2l [3J HI [5] because, even 
though their interactions are very weak, sizable amount of them may be produced in the 
early universe. If they are stable, they survive until today, resulting in the overclosure of the 
universe. Even if they are unstable, their lifetimes are often so long that they decay at a very 
late stage of the evolution of the universe. If they decay after the big-bang nucleosynthesis 
(BBN) starts, hadro- and photo-dissociation processes of light elements induced by their 
decay products may spoil the success of the BBN scenario. In order to avoid such a problem, 
the lifetime of X is required to be shorter than ~ 1 sec; otherwise, the abundance of X is 
bounded above [6j [7] . 

Even with lifetime shorter than ~ 1 sec, X may cause other cosmological difficulties. In 
our study, we consider one of such problems: the overproduction of the lightest superparticle 
(LSP) due to its decay. In many cases, X decays into superparticles which cascade down 
to the LSP. Then, if all the produced LSPs survive until today, it is often the case that the 
resultant mass density of the LSP becomes much larger than the present energy density of 
dark matter. In particular, the universe may be once dominated by X before the decay of 
X; in such a case, the overproduction of the LSP is a serious problem. 

If the pair annihilation cross section of the LSP is large enough, the LSPs produced by 
the decay of X may annihilate and the abundance of the LSP can be suppressed. In SUSY 
model, there exists a viable and well-motivated candidate of such LSP, which is neutral Wino 
W°. In models with the naive grand-unified-theory (GUT) relation among gaugino masses, 
Wino cannot be the LSP. However, even if the unification of the gauge group is realized, the 
GUT relation can be easily violated in various cases, like anomaly- mediation model [H [9], 
product-group unification [10], and so on. In particular, in the light of the recent discovery 
of the Higgs-like particle at the LHC [HI [12], the anomaly- mediation model with large scalar 
masses [8] (as large as ~ 10 — 100 TeV) is well- motivated because the relatively large SUSY 
Higgs mass of ~ 126 GeV can be realized if the scalar tops are so heavy [13] . It is well-known 
that Wino can naturally be the LSP in anomaly-mediation model. 

In this paper, we reconsider the possibility of the non-thermally produced Wino being 
dark matter. In particular, we precisely calculate the relic abundance of the neutral Wino. 
The possibility of Wino dark matter from the decay of exotic particles was first discussed in 
[SJ Q3] , in which it was shown that the relic density of Wino can be consistent with the dark 
matter density. Then, such a scenario has been applied to various cases rT5[ EH CEO ESI EES] • 
(For the case where the LSP is not Wino-like, see also [201 [211 [22j [23].) Compared to the 
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previous studies, we have carefully taken into account the following in the calculation of the 
relic Wino abundance: 

(i) Coannihilation effect among neutral and charged Winos, which becomes important 
when the decay temperature of X is higher than the mass difference between charged 
and neutral Winos. 

(ii) Sommerfeld enhancement of the pair annihilation cross section, which is large when 
Wino is heavier than ~ 1 TeV. 

The organization of this paper is as follows. In Section |2j we show relevant formulas to 
calculate the relic abundance of neutral Wino. In Section |3l we numerically calculate the 
relic abundance of Wino for several situations and discuss implications. Section @] is devoted 
to conclusions and discussion. 

2 Formulas 

In this section, we summarize the basic formulas to calculate the thermal relic abundance of 
the neutral Wino through the decay of a heavy particle X. 

Wino is SU (2) L-triplet, and there exist neutral and charged Winos, denoted as W° and 
W ± , respectively. In the present study, we consider the case where the neutral Wino is the 
LSP; in the following, we assume that superparticles other than Winos are so heavy that 
they are irrelevant at the time of the freeze-out of Winos. This is the case in the anomaly- 
mediated model, which is one of the important motivations of our study. In addition, because 
we assume that the Higgsinos are much heavier than Winos, the mass difference between W° 
and W ± are dominantly from one- loop diagrams with gauge bosons (i.e., 7, Z, and W ) and 
Winos inside the loop. Then, if the Wino mass is much smaller than the Higgsino mass, the 
mass difference is typically 150 — 170 MeV [2U [25], which is insensitive to the Wino mass. 

Because the mass difference is much smaller than the Wino mass, the number density 
of the charged Wino at the time of the freeze-out of Winos may be sizable. Thus, in the 
calculation of the relic density of W°, we include effects of all the possible channels of Wino 
annihilation. We consider the case where Winos are much lighter than other superparticles, 
so the coannihilation with those are neglected in the calculation of the relic density. 

Now we discuss the evolution of the number density of Wino. As we mentioned, we 
consider the case with a long-lived particle X which decays into MSSM particles. Using the 
fact that all the MSSM superparticles eventually decay into charged or neutral Wino, the 
relevant set of Boltzmann equations to calculate the relic abundance of Wino LSP is given 
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by 



din ~ 

-^L + 3Hn w = -(a cS v)(nl-n 2 Wcq ) + N w T x n x , (2.1) 
drix 



dt 



+ 3Hn x = -V x n X} (2.2) 



dt V 3 91nTy v r " u V 4<91nT 

where is the sum of the number densities of neutral and charged Winos, n x is the number 
density of X, and q is a heat injection into radiation as 

q = (m x - N w m w )T x n x +m w {a eS v)n 2 w . (2.4) 

In addition, p rad is the energy density of the relativistic component, and is related to the 
cosmic temperature T as 

PR = ^g*(T)T\ (2.5) 



where g*(T) is the effective number of relativistic degrees of freedomF 1 ! In our calculation, we 
approximated that the full particle content at the temperature above the QCD scale (which 
is taken to be 200 MeV in our analysis) is that of the MSSM, while that at the temperature 
below the QCD scale consists of photon, three generations of leptons, and pions. In addition, 
n w eq denotes the thermal-equilibrium value of n^, H is the expansion rate of the universe, 
Nyy is the averaged number of SUSY particles produced by the decay of one X, T x is 
the decay rate of X, and m x and are the masses of X and Wino, respectively!^! In 



the above Boltzmann equations, the thermally-averaged effective annihilation cross section 
(deffu) accounts both for the coannihilation effect and the Sommerfeld effect, which were not 
fully taken into accounts in previous analyses. 

The coannihilation processes are included by summing up the cross sections of all the 
relevant modes with appropriate weights: 

(a cS v) = ^riTjiaijv), (2.6) 



^ tl We use the fact that g* s (T) is numerically very close to g*(T), and approximate g* s (T) ~ g*(T) in our 
calculation, where g* s (T) is the effective number of masslcss degrees of freedom for the calculation of entropy 
density, which is related to the entropy density as 

2tt 2 , s „ 
s(T) = —g* s (T)T 3 . 

^ 2 Because we are interested in the case where charged and neutral Winos are quite degenerate, we denote 
the Wino masses as as far as we discuss the quantities which are insensitive to the mass difference. 
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where — W°, W + and W , and 

m 



n w 



(2.7) 



with rii being the number density of i. We assume that Winos are in kinetic equilibrium; 
when T <C m^, 



1 e -Am^/T 
M " = 1 + 2e -A%/T' ^ = ^ = 1 + 2e -Am^/T- ( 2 ' 8 ) 



If the mass difference between charged and neutral Winos Am^r is much larger than the 
background temperature T, only the lightest neutral Wino is relevant in the annihilation 
process. On the other hand, the coannihilation becomes effective for the temperature T > 
Am, 



For each annihilation process, thermally averaged cross section is obtained by# 3 

M = (£f f 2 / fvMe-™**'*. (2.9) 

Once the temperature of the universe decreases and the neutral and charged Winos become 
non-relativistic, the wave functions of annihilating Wino pairs are significantly deformed 
by the electroweak potential generated by the electroweak gauge boson exchanges. The 
resultant annihilation cross sections, (TijV, are significantly enhanced or suppressed due to 
the Sommerfeld effect [28| [29| 130] . The Sommerfeld effect is more important for larger Wino 
mass, > 1 TeV, because the electroweak potential behaves as a long range force in such a 
mass region. As a result, the thermally-averaged cross sections show non-trivial dependence 
on the cosmic temperature T and the Wino mass m^. Two-body system of Winos can be 
classified by the quantum numbers Q (electric charge) and S (spin), and the Sommerfeld 
enhancement factors are evaluated for fixed values of these quantum numbers. For each set 
of (Q,S), possible decay modes and the decay widths are summarized in Appendix. (See 
Tabled!) 

In Fig. 121 we show the effective annihilation cross section for several cosmic temperatures, 
T = 10 _5 7n^, 10~ 4 m l y, lO^m^, and T 3> m^. (Here and hereafter, we take Am^ = 
170 MeV in our numerical calculation.) For comparison, we also show the annihilation cross 
section of neutral Wino without the Sommerfeld effect in the same figure. Obviously, the 
cross section is significantly enhanced in particular for the Wino mass larger than ~ 1 TeV, 
and the enhancement factor becomes larger for lower temperature. In addition, we can see 
the resonance structure at around ~ 2.4 TeV, which occurs due to the existence of zero- 
energy bound states. The precise position of the resonance is determined by the structure 



# 3 The energetic Winos produced by the X decay show non-trivial velocity distribution at first. However, 
charged Winos soon lose their energies through the electromagnetic interactions with the thermal background, 
and also sizable fractions of neutral Winos can be thermalized through the inelastic interactions [27] . 
Then, here we simply assume the Boltzmann distributions to obtain the thermally-averaged cross section. 
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Figure 1: Thermally averaged total cross section (a e gv ) for T = lO^^m^ (pink line), lO^m^ 
(blue line), I0~ 3 myy (green line), and T 3> (red line). For comparison, we also show the 
annihilation cross section of neutral Wino without the Sommerfeld effect (black line). The 
vertical axis is the the Wino mass, and we take Am^ = 170 MeV. 

of electroweak potential, and depends on the mass difference between charged and neutral 
Winos, Arriyy. For larger Am^, the resonance peak is shifted to the heavier Wino mass. 

As is seen from the Fig. [2j the annihilation cross section at the decay temperature can be 
significantly enhanced by the Sommerfeld effect. We stress here again that once coannihila- 
tion and Sommerfeld effects are included, the cross sections show non-trivial dependence on 
T. In order to precisely take into account these effects, we solve the Boltzmann equations 
numerically. 

The relic density of the neutral Wino can be calculated by solving Eqs. ( 12. II) — ( 12. 31) with 
relevant initial condition. The initial values of n^, rix, and p ra d depend on cosmological 
scenarios and the properties of X. In the next section, we consider several well-motivated 
scenarios and calculate the relic density by numerically solving the Boltzmann equations. 

3 Relic Abundance of Neutral Wino 

3.1 Case with X domination in the early universe 

First, we consider the case where the particle X once dominates the universe. In such a case, 
at the cosmic time t <C T^ 1 , the energy density of X is much larger than that of radiation, 
while the energy density of radiation (i.e., so-called "dilute plasma") scales as a -3 / 8 (with a 
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being the scale factor). With such an initial condition, we solve the Boltzmann equations 
to calculate the relic abundance of W°. We assume that the maximal temperature of the 
dilute plasma is much higher than the Wino mass; then, the pair annihilation rate is initially 
much larger than the expansion rate of the universe. In such a case, the production and 
annihilation terms in the right-hand side of Eq. (12.1 ft (almost) balance, so we take the initial 
value of the number density of Wino as n|™ lt ' ) = ^Jn 2 ^ ^ + N-^T x^x^^v)" 1 . (However, the 

resultant relic abundance is insensitive to the initial value of as far as the initial condition 
is set at the cosmic time with the background temperature much higher than m^.) With the 
initial condition given above, we calculate the number density of Wino after the completion 
of the decay of X, and evaluate the yield variable of Wino 

Y W = ^, (3.10) 

s 

where s is the entropy density. Using the fact that becomes constant of time at low 
enough temperature, we calculate the density parameter as 

-l 

Pcrit 



n w = m w [Y w ] t>>r -i , (3.11) 

x V-Snow/ 

where p crit is the critical density of the universe while s now is the present entropy density, and 
their ratio is given by p cr it/s now ~ 3.6/i 2 x 10~ 9 GeV, with h being the Hubble constant in 
units of 100 km/sec/Mpc. (In our numerical calculation, we use h = 0.697 |31j.) Comparing 
with the present density parameter of dark matter, we derive constraints on the model 
parameters. We use the following value of the dark matter density as the canonical value 

Q c h 2 = 0.1146. (3.12) 

In the case where X once dominates the universe, the relic density of Wino depends on 
the following parameters: T x (decay rate of X), m x (mass of X), (Wino mass), and 
Nyy (averaged number of Winos produced by the decay of one X). In particular, the relic 
density depends on the decay rate Tx- To discuss the dependence on Tx, it is convenient to 
define the "decay temperature" as 



10 



1/4 



T x = ———Mife , (3.13) 



g^TxW 

where Mpi ~ 2.4 x 10 18 GeV is the reduced Planck scale. Notice that Tx corresponds to the 
cosmic temperature at the time of X decay. 

In Fig. |2l we show the contours of constant Sl^h 2 for the case with = 1. (In this 
subsection, we adopt the relation mx = 5to^; the resultant value of VLy^h 2 is insensitive to 
this relation as far as the pair annihilation of Winos once becomes effective.) As one can see, 
flyyh 2 is insensitive to T x when the decay temperature is high, while fly^h 2 increases as T x 
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Figure 2: Contours of constant fl^h 2 for = 1 and m x = 5m^ for the case where X 
once dominates the universe. Numbers in the figure are Q^h 2 . 



becomes lower when the decay temperature is relatively low. Such behaviors can be easily 
understood. With high enough decay temperature, Winos are still in chemical equilibrium 
after the completion of the X decay, and hence the relic abundance of Wino is given by the 
thermal relic density. On the contrary, with low decay temperature, Winos cannot be in 
chemical equilibrium. Even in such a case, the pair annihilation of Wino proceeds as far as 
the annihilation rate is larger than the expansion rate of the universe. In such a case, the 
yield value of Wino is approximated by 



Y w ~ min 



3T X N w n x (T x ) 



{a eS v)s(T x ) s(T 



Xj 



(3.14) 



where n x (T x ) represents the number density of X just before the decay. In the present case, 
the first term, which is approximately proportional to T x , is smaller than the second one, 
resulting in the enhancement of fly^h 2 with lower decay temperature. 

In order to realize Wino dark matter scenario (i.e., = f) c )? relatively high value of 
T x is needed. For m w = 300 GeV (500 GeV, 1 TeV, 2 TeV), T x = 1.3 GeV (4.0 GeV, 
23 GeV, 150 GeV) is necessary. If T x is lower, the universe is overdosed. This has an 
important implication when the cosmological moduli fields play the role of X, as we discuss 
below. In addition, we can see a significant suppression of Sl^h 2 when ~ 2.4 TeV. Such 
a suppression is due to the significant enhancement of the annihilation cross section by the 
Sommerfeld effect. 
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Figure 3: ^ly^h 2 as a function of for = 1 and T x = 10 GeV with Sommerfeld 
and coannihilation effects (red line), with coannihilation effect only (green line) and without 
neither Sommerfeld nor coannihilation effects (blue line). 

In order to see the importance of Sommerfeld enhancement, we also calculate fl^yh 2 
neglecting the Sommerfeld enhancement. In addition, for comparison, we also show the 
result without the effect of coannihilation, for which we use (a e sv) = (<Twowo^w+w- v ) 
(without Sommerfeld effect). The results are shown in Fig. [3] for Nyy = 1 and Tx = 10 GeV. 
The Sommerfeld effect significantly changes fl^h 2 when ^ 1 TeV. 

In the case where X dominates the universe, the N-^ dependence of VL-^h 2 is very weak 
unless is significantly suppressed. In Fig. 0], we plot contours of constant fly^h 2 with 
Nyy = 10~ 6 . Even with such a small value of N^, the contour of Qyy = Q c does not change 
much compared to that given in Fig. [2j Thus, unless <C 10 -6 , Tx should be higher than 
0.6 GeV (2 GeV, 10 GeV, and 90 GeV) for m w = 300 GeV (m w = 500 GeV, m w = 1 TeV, 
and m^y = 2 TeV). In Fig. EJ we show the contours of the upper bound on Nyy on the 
vs. T R plane, requiring < Q c . As a result, we can see that < O(10 -8 ) is required if 
Tx is so small that the pair annihilation of Winos is ineffective. 

One of the important applications of our study is to the case with weakly interacting 
late decaying scalar fields, like scalar fields responsible for the SUSY breaking (i.e., so- 
called the Polonyi field) and the moduli fields in string theory. (Hereafter, we call them 
"moduli fields.") Moduli fields are expected to have interactions suppressed by the Planck 
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Figure 4: Same as Fig. [2} except for = 1CT 6 . 



scale. A modulus field may have a large initial amplitude. Because of the weakness of the 
interaction, the lifetime of the modulus may be so long that its decay occurs at a late stas^e 
of the evolution of the universe. It is well known that the late-time decay of the modulus 
field is dangerous because it destroys the light elements produced by the BBN reactions, 
resulting in the spoil of the success of the BBN scenario. Such a problem can be avoided 
if the decay rate of the modulus is somehow enhanced so that the lifetime of the modulus 
becomes shorter than ~ 1 sec. Even so, the decay of the modulus field produces significant 
amount of superparticles, which may result in the overproduction of the LSP. With Wino 
LSP, this problem may be avoided. 

The decay rate of a modulus field depends on how it interacts with the MSSM fields. For 
example, it may interact as 

Ant = / d 2 6XW a W a + h.c, (3.15) 



Mi 



pi 



where the "hat" stands for superfields, VV is the gauge field strength superfield, and \q is a 
coupling constant. Then, the scalar component X decays into a pair of vector boson V with 
the following decay rate: 



NjXq m 3 x 

x^vv = — 7T9~? (3.16) 



4tt M|j : 



where N{ is the number of the possible final states; for example, N{ = iV 2 — 1 for an SU(iV) 
gauge group. (Here, we assume that the vector boson is much lighter than X.) In addition, 
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Figure 5: Contours of Qyy = Q c for = 1 (blue line), Nyy = 10~ 6 (green line) and 
N w = KT 8 (red line). 

with the interaction given in Eq. A3. 1 5j) . X may decay into a pair of gauginos A. The 
decay rate is model-dependent, and is typically Tx^aa ~ ^x^vv [321 El]- Approximating 
Tx ~ Tx-^yy, the decay temperature is estimated to be 

T x ~ 0.01 GeV x A G ( mx Y * , (3.17) 
G V100TeV/ ' V ' 

where we have used Nf = 12, which is the number of the gauge bosons in the standard 
model. For = 300 GeV (1 TeV), the Ac-parameter is required to be larger than 430, 
130, and 42 (7600, 2300, and 74) for m x = 10 TeV, 100 TeV, and 1000 TeV, respectively, 
where we have assumed Nyy ~ 1. It seems that if such a modulus with Ac ~ 0(1) once 
dominates the universe it tends to cause overproduction of Winos. In order for Winos not 
to overdose the universe, large mx or enhanced Aq is needed. 



3.2 Case with gravitino 

Next, we consider the case where gravitino plays the role of X. Gravitino may not dominate 
the universe, but a significant amount of W may be produced by its decay. 

Gravitino is the superpartner of graviton, and it couples to the supercurrent so it interacts 
with all the supermultiplets. If unstable, gravitino decays into an ordinary (visible-sector) 
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particle and its superpartner. Thus, decay of gravitino results in the production of the 
LSP because all the produced superparticles cascade down to the LSP. The interactions of 
gravitino are suppressed by inverse powers of the Planck scale so that the decay rate of 
gravitino is extremely small. It is well known that the late-time decay of the gravitino may 
spoil the success of the BBN if the lifetimes of the gravitino is longer than ~ 1 sec [HJ 171 132]. 
In order to avoid such a problem, we concentrate on the case where the lifetime is shorter 
than ~ 1 sec, which is realized if the gravitino mass is larger than ~ O(10 TeV). Assuming 
that the gaugino masses are at the TeV scale and are much smaller than the gravitino mass, 
the lifetime of the gravitino is estimated as 



T3/2 



0.4 sec x N^ 1 r ^,, 1 3 , (3.18) 
G V 100 TeV/ v ' 



where 7/13/2 is gravitino mass. In addition, Nq is the number of gauge multiplets to which 
gravitino decays. (In the following numerical study, we take the MSSM value of Nq = 12.) 
We also assume that the superparticles other than gauginos are as heavy as gravitino so that 
the decay modes into those superparticles are negligible. 

Even though gravitino is very weakly interacting, gravitinos are produced by scattering 
processes of particles in thermal bath. The abundance of gravitino (before its decay) is 
approximately proportional to the reheating temperature after inflation; for the case where 
the gravitino mass is significantly larger than gaugino masses, the yield variable of gravitino, 
which is defined as Y 3 / 2 = n 3 / 2 /s with n 3 / 2 being the number density of gravitino, is given 

by mB 



[V 3/2 ] ^ 2.3 x 10~ 14 x T^ 8) 1 + 0.015 In - 0.0009 In 2 T$ } , (3.19) 



I t<r 3/2 



(8) 



where = Tr/10 8 GeV, with Tr being the reheating temperature after inflationj^j which 
is related to the decay rate of the inflaton r inf as 



10 ^ 1/4 



T ^{^P M " Tl ') ■ (3 ' 20) 

We calculate the density parameter of Wino by numerically solving the Boltzmann equa- 
tions given in Section [2j Because the primordial abundance of the gravitino is approximately 
proportional to the reheating temperature, fl^ also increases with higher reheating temper- 
ature. 

In Fig. [6j we show the contours of fl^ = fl c for several values of Tr on vs. m 3 / 2 
plane 1^1 The contour of Vt^ = Vt c is almost independent of the reheating temperature if 



# 4 If there exists condensation of an exotic scalar (like inflaton, moduli fields, and so on), gravitino may 
be produced by the decay process [33l [Ml ES] ■ We assume that the abundance of the gravitino from such a 
decay process, which is highly model-dependent, is negligible. 

# 5 Here, we assume that there is no significant entropy production after inflation. 

# 6 In the anomaly-mediation scenario with ~ O(100 GeV — 1 TeV), the gravitino mass is expected to 
be of O(100 TeV). Here, however, we vary 7713/2 up to O(10 5 TeV) to provide information about the case 
with extremely heavy gravitino. 
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Figure 6: Contours of Vt^ = Q c for several values of the reheating temperature. Here we 
take Nq = 12. Numbers in the figure are Tr in units of GeV. 



Tr is high enough. This is due to the fact that, with high enough reheating temperature 
and large enough gravitino mass, the Wino abundance at the time of the gravitino decay 
is so large that the pair annihilation of Wino becomes effective. Then, as discussed in the 
previous subsection, Qyy is determined by the decay temperature which is determined by 
7713/2 and becomes insensitive to the primordial abundance of the gravitino. On the contrary, 
if the gravitino mass is smaller than ~ 10 3 TeV, the contour of Qyy = Q c becomes insensitive 
to the gravitino mass. This is because, if the gravitino mass is small, the decay temperature 
of the gravitino becomes so low that the annihilation of Wino is inefficient. Then, almost all 
Winos produced by the gravitino decay survive until today, and the relic density of Wino 
is approximately given by the sum of thermal relic density and the non-thermal one from 
the gravitino decay: Y^ ~ Y^ + p3/2]t<T 3/2 , with being the thermal abundance of 
Wino. The thermal relic density becomes larger than the present dark matter density if 
m w ~ 2.9 TeV, so the Wino mass larger than ~ 2.9 TeV is forbidden irrespective of Tr. On 
the contrary, with smaller Wino mass, we obtain upper bound on the reheating temperature 
by requiring Q-^ < Q c , using the fact that the primordial number density of gravitino is 
approximately proportional to Tr. Such a bound is shown in Fig. [7J for several values of 
gravitino mass. 

One implication of our result is on the leptogenesis scenario [36J. The leptogenesis sce- 
nario requires the reheating temperature to be higher than ~ 10 9 GeV in order to generate 
large enough amount of the baryon asymmetry of the universe [371 138] . We have seen that 
such a high reheating temperature can be realized in the Wino LSP case. 
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Figure 7: Upper bound on the reheating temperature as a function of Wino mass for several 
values of the gravitino mass. Numbers in the figure are gravitino mass in units of TeV. (The 
bound for m 3 / 2 = 10 2 TeV (blue line) almost overlaps with that of m.3/2 = 10 TeV.) 



3.3 Case with SUSY Peccei-Quinn model 

The Peccei-Quinn (PQ) mechanism [39l HO] is an attractive solution to the strong CP prob- 
lem. If we embed the PQ mechanism into supersymmetric model, there exists the fermionic 
superpartner of the axion called axino. Axinos are copiously produced in the thermal bath 
of early universe and decay into MSSM particles at a late stage of the cosmic expansion. 
In this subsection, we consider the case where the axino a plays the role of X, which is 
embedded in the axion multiplet A aa^ 7 

A= - r (a + ia) + V29a + (F-tevm). (3.21) 
v2 

(In this subsection, the field X is denoted as a.) 

To make our discussion concrete, we assume that the PQ charges of all the MSSM fields 
are zero, and that the PQ fermions are embedded into full multiplets of SU (5) grand unified 
gauge group. In such a model, denoting the PQ scale as f a , the axion multiplet A has the 

# 7 The real scalar component of A (i.e., saxion) may also cause cosmological difficulties. Cosmology with 
saxion depends on the model of PQ sector as well as on thermal history. See [4TJ [42J H31 EJ [45] for more 
detailed discussion. 
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following interaction term: 

+ h.c., (3.22) 



4V2vr/ a 



a 3 AG a g a + a 2 AVV a VV a + \a x AE a E 



where Q, VV, and B are field strength superfields for SU(3)c, SU(2)i, and C/(l)y, respec- 
tively, and a a = g 2 /47r with g a (a = 1 — 3) being gauge coupling constants. (The summation 
over the gauge indices are implicit for SU(3)c and SU{2) L .) 

As in the case of gravitino, the primordial abundance of the axino depends on the reheat- 
ing temperature after inflation [25]. We adopt the yield variable of axino (i.e., Y"„ = n„/s, 
with ria being the number density of axino) evaluated in |47j : 



[^Wrr 1 - min 



F 5 (eq) ,0.20 x a 3 In 



0.0977 \ / Tr \ / f a V 

y Vio 7 Gevy Vio 11 Gevy 



«3 



(3.23) 



where ~ 1.8 x 1CT 3 is the thermal abundance of axino. If the axino is stable, the 

axinos produced in the early universe survive until today. Then, in order not to overdose 
the universe with the mass density of axino, severe upper bound on the reheating temperature 
is obtained [I5J H?l HH H9]. Even if all the LSPs produced by the axino decay remain until 
today, which is the case if the annihilation cross section of the LSP is small, a stringent 
upper bound on Tr still exists. With a large annihilation cross section of dark matter, this 
problem may be avoided [50] . 

One of the important parameters to calculate the relic Wino abundance produced by the 
axino decay is the decay rate of axino. In the present setup, the axino decays dominantly 
into gauge boson and gaugino pair; with the interaction terms given in Eq. (13.221) . the decay 
rate is given by 

Ta = Ta^gg + V^ w ±w^ + I 'g^jyo + ^5,-¥jW° + ^a-*ZB + ^a-^Bi (3.24) 

where, if kinematically allowed [IT] . 

8a 3 2 ma 3 , 2 x 3 , Q _s 

T~ a ^ w±WT = 128 ^-3 ^"2 K (ywiVw), (3.26) 



Op 1 COS^ 0\zy TTt q " 

T ~a^zw° = 128tt 3 f a 2 K (yw>yz), (3-27) 



a 2 2 sm 2 9 w m a 3 , 2]3 
l a-^ 1 W° ~ 128vr 3 ~Y T " w > ' ^ ' 



5\ 2 afsin 2 ^ m a 3 

3j 128tT 3 X 

2 „,2_-.„2/l ™ 3 



r~ri = (?) (3-29) 



/5\ o^cos 2 ^ m a 3 , 2 \3 / q qn v 

r ^ = UJ "t^t^ ^ } ' (3 - 30) 
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where 6w is the Weinberg angle, yj = mj/rria denotes the mass of the particle I normalized 
by the axino mass m a , and 



K(Vi, V2) = a/1 + Vx + Vi ~ 1y\ - 22/1 " Hyl ' 



.(3.31) 



The decay temperature, and hence the resultant Wino abundance, depend on the mass 
spectrum of superparticles, in particular, that of gauginos. In the case where the gluino is 
lighter than the axino, the axino dominantly decays into gluon and gluino pair. In such a 
case, the decay rate of the axino becomes relatively large. In the opposite case, the dominant 
decay modes of the axino are a — > W ± W Zf , a — > ZW°, and a — > jW°, and the decay rate of 
the axino is suppressed. So, in the latter case, VL^ becomes larger compared to the former 
case. 

In Fig. |HJ we show the contours of constant VL^h 2 for the case with high enough reheating 
temperature ([Y" s ], „ r -i = ), taking : : m„ : m a = 1 : 3 : 7 : 10. (The gaugino 
masses are assumed to obey the anomaly- mediation relation J8j [9].) With such a mass 
spectrum, the axino dominantly decays into gluon and gluino pair. Then, Qyy = Q c is 
realized with f a = 1.9 x 10 11 GeV (1.3 x 10 n GeV, 6.3 x 10 10 GeV, 2.7 x 10 10 GeV) for 
m w = 300 GeV (m w = 500 GeV, m w = 1 TeV, and m w = 2 TeV). 

If we consider the case where the decay mode into the gluon and gluino pair is kinemat- 
ically blocked, the decay rate is suppressed. In such a case, the PQ scale which realizes the 
Wino dark matter becomes smaller. In Fig. [91 we show the result taking : : m~ g : 
rria = 1:3:7:5. Then, the value of f a giving rise to fl^ = Q c is given by 1.5 x 10 10 GeV 
(1.0 x 10 10 GeV, 5.1 x 10 9 GeV, 2.2 x 10 9 GeV) for m w = 300 GeV (m w = 500 GeV, 
m^y = 1 TeV, and = 2 TeV). 

We have also studied how the required value of f a to realize Wino dark matter depends 
on the axino mass. In Fig. [lOj we show contours of constant f a which gives Qyy = Q c for 

[^Wr- 1 = Y t q) on m w vs - m ~al m w plane. 

a 

So far, we have adopted the thermal abundance of axino. Even if the primordial abun- 
dance of axino is smaller, the resultant Wino density does not change as far as Winos 
produced by decay of axino are much enough to pair-annihilate. We also performed the 
calculation with lower reheating temperature, and checked that the results are more or less 
unchanged. 



4 Conclusions and Discussion 

In this paper, we have studied the possibility of the Wino cold dark matter in supersymmetric 
models, paying particular attentions to the scenario in which the decay of a long-lived exotic 
particle X produces significant amount of Wino LSP. We have numerically calculated the 
relic abundance of the neutral Wino, carefully taking account of the effects of coannihilation 
and Sommerfeld effects. We have seen that the Sommerfeld effect drastically enhances the 
pair annihilation cross section of Wino if > 1 TeV. We have studied the cases where 
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Figure 8: Contour of fly^h 2 of the present universe for mass ratio : : m g : m a = 1 
3 : 7 : 10. Here, we take [5 / a] t< rr 1 = ■ Numbers in the figure are the values of ft^h 2 . 
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Figure 9: Same as Fig. El except for rriy, : : m g : m a = 1:3:7:5. 
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Figure 10: Contours of constant f a which gives Vt^ = tt c for \Y & ]. <r -i = 3^ . The 
horizontal axis is the Wino mass while the vertical one is the ratio m^/m^. The gaugino 
masses are assumed to obey the anomaly-mediation relation :m^: rrig = 1:3:7. 



moduli fields, gravitino, or axino plays the role of X, and derived the constraints on the 
model parameters in each case. 

So far, we have not discussed cosmological and astrophysical constraints on the Wino 
dark matter. Because the neutral Wino has large annihilation cross section, the Wino dark 
matter scenario conflicts with astrophysical and cosmological constraints if the Wino mass is 
too small. One of the constraints is from the negative observation of high energy 7-ray from 
dwarf galaxies, from which the mass regions < 400 GeV are disfavored (511 [52] although 
astrophysical uncertainties exist. Another constraint is from BBN. The relic Wino may pair 
annihilated during and after the BBN epoch, which causes photo- and hadro-dissociation 
processes. In order not to spoil the success of the BBN scenario, the Wino mass is required 
to be larger than ~ 200 GeV [53J. Thus, these astrophysical and cosmological constraints 
do not exclude most of the parameter regions we have studied. 

Acknowledgements: This work is supported in part by Grant-in-Aid for Scientific research 
from the Ministry of Education, Science, Sports, and Culture (MEXT), Japan, No. 22244021, 
No. 23104008, and No. 60322997. 
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A Sommerfeld Enhancement Factors 



Here we present the possible decay modes and the decay widths for each set of (Q, S), which 
are summarized in Table [U The annihilation cross section to the / /' final state is obtained 
by 



av = c- 



Ti 2 a 2 2 



m 



^(2S + l)P(x f ,x f )D\A(< 



(A.1) 



W 5=0,1 



where c = 2 for the annihilation of identical particles, otherwise c = 1. 
The phase factor P(xi,x 2 ) is given by 



/ / y» 2 | ry>2 ( ™ 2 ,-y»2\2 

P( Xl ,x 2 ) = sjl - + {Xl 16 X2} . 



In addition, the functions fi(x) are defined as 

4(1 - 40(4 + 20x 2 w + 340(2 + x 



fi(x w , x z ) 

f2(x W ,X Z ) 

fs(,x,x h ) 
U(x z ,x f ) 

f 6 (x w ,x z ) 



2 

w 



X 



2 \2 



(2-x^(4-4)2 
{(4 -x 2 w - x 2 z f - Ax 2 w x%} (16 + 40x 2 w + 40x 2 z + 10x 2 w x 2 z + x 4 w + x\ 



(4-x 2 w ) 2 (4-x 



1 - 



1 - 



X 

7 



2 \ - 2 

xP 



1 + 



P\ ry 2 ,-yi 2 



+ 
,,2 



(x 2 - 4) 



w 

2\2 



X 



2 \2 



16 



2 / 4 

1 ~~ c / x z + ~2" x z 



^ (1 + 2 C/ 4 - cjx|) 



2 \ -2 



,7-> 2 I ™ 2 



32 



2 1 2 \ -2 



rf> 2 I ,-y»2 

•^vy ^ -^z 



+ 



.7; 



X 



2 \2 
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where c/ = <5/Svf/(2T 3 ^ l ) for the particle with the electric charge Q/ and the weak isospin 



J 3L- 
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Table 1: Table of parameters to calculate annihilation cross sections for each initial and final 
states. Tree-level cross sections are recovered by setting Af^ ^ = A^ ^ = A^ 0,1 ^ = A^ 1,0 ^ = 
,4(1,1) = ,4(2,0) _ 1 anc [ A21 = Af^ = 0. Here, xj is the mass of the particle / normalized 
by rriyy, xi = Tni/m^. f denotes quarks and leptons. Nc is the color factor; Nc = 3 (1) for 
quarks (leptons). For the calculation of the Sommerfeld factors A^' s ^ for states with charge 
Q and spin 5*, see [30J. 
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